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SEMI-REGULAR TESSELLATION TYPE GRAPHS 



EUGEN J. IONASCU 

Abstract. We study the problem of half-domination sets of vertices in transitive infinite graphs 
generated by regular or semi-regular tessellations of the plane. In some cases, the results obtained 
are sharp and in the rest, we show upper bounds for the average densities of vertices in half- 
domination sets. 



1. Introduction 

By a tiling of the plane one understands a countable union of closed sets (called tiles) whose 
union is the whole plane and with the property that every two of these sets have disjoint interiors. 
The term tessellation is a more modern one that is used mostly for special tilings. We are going 
to be interested in the tilings in which the closed sets are either all copies of one single regular 
convex polygon (regular tessellations) or several ones (semi-regular tessellations) and in which each 
vertex has the same vertex arrangement (the number and order of regular polygons meeting at a 
vertex). Also we are considering the edge-to-edge restriction, meaning that every two tiles either do 
not intersect or intersect along a common edge, or at a common vertex. According to [I], there are 
three regular edge-to-edge tessellations and eight semi-regular tessellations (see [3], pages 58-59). 
The generic tiles in a regular tessellation, or in a semi-regular one, are usually called prototiles. 
For instance, in a regular tessellation the prototiles are squares, equilateral triangles or regular 
hexagons. We will refer to these tessellations by an abbreviation that stands for the ordered tuple 
of positive integers that give the so called vertex arrangement (i.e. the number of sides of the regular 
polygons around a vertex starting with the smallest size and taking into account counterclockwise 
order). The abbreviation is usually using the convention with the powers similar to that used in 
the standard prime factorization of natural numbers. So, the regular tessellation with squares is 
referred to as (4 4 ). We refer to Figure Q] for the rest of the notation. 

Each such tessellation is periodic in the sense that there exists a cluster of tiles formed by regular 
polygons, which by translations generated by only two vectors, say v\ and t>2, covers the whole 
plane, and the resulting tiling is the given tessellation (Figure [1]). 
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(3 4 ,6) 

Figure 2. GW 6 w 3 



For each such tessellation T, we associate an infinite graph, G-7-,00 in the following way. For 
each regular polygon in the tessellation we have a vertex in G-j-,oo, and the edges in this graph are 
determined by every two polygons that share a common edge within the tessellation. This way, for 
instance, for the regular tessellation with squares (4 4 ), we obtain what is usually called the infinite 
grid graph. Also, if m, n £ N we can construct a graph Gf,mn in the same way as before using 
only m copies of the cluster of tiles generating T shown in Figure [T] in the direction of v\ and n 
copies in the direction of v%. For T = (3 4 , 6), m = 4 and n = 3 we obtain the graph generated by 
the polygons in the portion of the tessellation T shown in Figure [2j 

Clearly, the choice of the tiles on Figured] is not unique but whatever one takes for these tiles it 
is not going to be relevant in our calculation of densities. For each such graph as before, and one 
of its vertices v , let us denote by d(v) the number of adjacent vertices to v, known as the degree of 
the vertex v. We define a set of vertices S to be half- dependent if for each vertex v G S the number 
of adjacent vertices to v that are in S is less than or equal to [~^J • 

Let m, n £ N be arbitrary, and for each graph GT,m,n we denote by PT,m,n the maximum 
cardinality of a half-dependent set in Gj^ m ^ n) divided by the number of vertices of G-j-,m,n- Hence 
one may consider the number 
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(1) PT = lim SUp PT,m,n, 

m,n— too 

which represents, heuristically speaking, the highest proportion in which one can distribute the 
vertices in a half-domination set. For instance we will show that in the case of a regular tessellation 
with regular hexagons (see Figure [3]) (b) the number defined above is |. 




(a) (b) 

Figure 3. Regular tessellations with squares and regular hexagons 



In the regular tessellation with squares we have shown in [7J that the number defined in (TjQ) is also 
| (Figure [3ja)). In this paper we are interested in the values of pf for these types of tessellations. 

2. Various techniques 

First we are going to use one of main ideas in [B] and [7], and some classical linear optimization 
techniques. 

2.1. Integer Linear Programming Systems. The regular tessellation (6 3 ) can be treated simply 
in the following way. We are going to work with the graph Cr6 3 ,n,n obtained from translating n- 
times a prototile so that the each two neighboring tiles share an edge and then translate the whole 
row of n-tiles in such a way that two neighboring rows fit together to give the tessellation of a 
rhombic n-by-n region. We assign to each tile a variable Xij which can be or 1: if a tile is part of 
the half-dominating set (colored blue in Figure [3b) of maximum cardinality then its variable x% j 
is equal to 1, and if it is not, Xij = 0. For most of the vertices of G 6 3 nn the degree is 6 so each 
vertex in the half-dominating set, denoted by V, must have at most 3 other vertices adjacent to it 
which are in V. Let us denote by N(i,j) the adjacent indices to the vertex indexed with (i,j). We 
can write the half-domination condition as 
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X tj := ^2 Xk > 1 ~ 3 for a11 ^'j) With 1 < ^ ^ < n ' (*>j) G 1/ - 

(k,i)eN{i,j) 

For a vertex not in F we simply have no restriction on the above sum (maximum is 6). The trick 
is to write an inequality that encompasses both situations. In this case, we see that the following 
inequality accomplishes exactly that: 

3x m + Xk > 1 — ® ^ or a ^ w ^h 1 < i,j < n. 

(k,l)EN(i,j) 

For the vertices on the boundary we have similar inequalities. Summing all these inequalities up, 
we obtain, 3|V'| + 6|V| + E < 6n 2 , where V are the vertices in V in the interior and S are the 
number of those on the boundary. Since clearly lim 4 = we see that lim^i < §. In Fi gure [3)3, 
we exemplify an arrangement which accomplishes this density. Hence, p G 3 = |. 

In the case of the tessellation with equilateral triangles, the inequality above changes into 

2xa + x/x' < 3 for all As not at the boundary. 

A'eJV(A) 

The argument above gives the inequality /9 3 e < |. The arrangement in Figured] has a density which 
is only j§. In this case it is not possible to achieve the density |. 

2.2. Toroidal type graphs. One way to avoid to deal with the boundary tiles, is to form toroidal 
type graphs obtained from tessellations. We have shown in [6] that this does not affect the maximum 
density, in the sense that both situations tend at the limit to the same density value. Let us begin 
with the case T = 3 6 . For n € N, we take the parallelogram which gives the graph G36 nn 
and identify the opposite edges (without changing the direction). This gives rise to a similar 
graph which we will denote by T 3 e n . This graph is regular and the equations we get can be easily 
described and implemented in LPSolve. Let us introduce the variables a^fc, i, j € {0, 1, 2, n — 1} 
and k G {1,2}, in the following way. The index refers to the translation of the minimal 

parallelogram (Figure [1]) i places in the horizontal direction and j places in the 60° direction. The 
index k corresponds to the lower (k = 1) and upper equilateral triangle within this parallelogram 
(k = 0). We get a number of 2n 2 vertices of this regular graph. The optimization equations are 
simply 
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(2) 



' 1x%^\ + Xi,j t 2 + x u,j,2 + Ki,v,2 < 3, where u = i — 1 (mod n), 

v = j — 1 (mod n), u, t> € {0, 1, 2, n — 1} 
% x i,j,2 + + x mj,i + x i,v,i < 3, where w = i + 1 (mod re), 



k i = j + 1 (mod n.),io,t G {0, 1,2, ...,n - l},z, j G {0, 1,2, ...,n - 1}. 
We denote by 036^ the best densities for this graph. Observe that p 3 e n < p^, m if n divides m. 
Also, as proved in [6], we have p 3 e = supp 3 6 n . We have the following densities: 
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16 


14 

25 


5 
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- 49 
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- 16 


> k 
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Density 

Figure 4. Regular tessellation with equilateral triangles 

Unfortunately, LPSolve takes too long to finish the analysis for n > 7 (about 100 variables). In order 
to reduce the number of variables, and of course taking advantage of the rotational symmetry, one 
can try to do a different matching of the boundaries as in Figure [5l This corresponds to rotations 
of 180° around the midpoints of the sides for the basic triangle of tiles ABC. For instance, the 
corner tiles are all adjacent to each other. 

The equations that we need to use, are those given by ([2]) and in addition 

x i j i = x uv 2 where u = n — 1 — i (mod n), and v = n — 1 — j (mod n), 

(3) 

u,v e {0, 1,2, ...,n - l},for all i,j £ {0, 1,2, ...,n - 1}. 
We will refer to the graph obtained by these identifications as K^e n . This graph is a regular graph 
with n 2 vertices and degree 3. All the densities are the same as in the table above but LPSolve 
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A=B 




Figure 5. Klein type toroidal identification 

finishes in a reasonable time and the arrangement in Figure U] satisfies this new restriction. In the 
case n = 4 we can see that the arrangement in Figure 0] is also the best because if we could add 
another triangle, that would increase the density to | > |, which we know is not possible by the 
upper bound established earlier. It is important that the best density in the case K^a 8 is the same 
as in the case 4 . We conjecture that the best density is given by such a matching of two toroidal 
graphs one having the dimension double the dimension of the other. 

2.3. Upper bounds. For the problem above we can adopt the method of weighted objective 
function as described in [TJ. We used the weights all equal to 1 in the interior and zero on the 
boundary for i^s 13 . The upper bound obtained is p 3 6 < w 0.59375. This is within « 
0.0160124 to jq which makes it plausible that p 3 6 = yg. 

For the graph K 3 6 n , we may add to the system of inequalities (j2J) and the conditions < 
x i,j,i, x i,j,2 < 1 and Y2i j{ x i,j,i + Xi,j,2) = k. These inequalities describe a polytope in n 2 dimensions. 
Finding the maximum cardinality of a half dependent set is equivalent to finding the smallest k for 
which there is no lattice point in the corresponding polytope. There exists a theory which counts 
the number of lattice points in polytopes which was started by Eugene Ehrhart in 1960's (see [1] 
and [2]). The theory simplifies significantly if the polytope has integer vertices. Unfortunately our 
polytope has rational vertices. Theoretically, there exists a quasi-polynomial P(t, k) of degree n 2 
which counts the lattice points contained in the dilation of the polytope by t. We want the smallest 
k n such that P(l,k n + 1) = 0. Let us make the observation that P(l,k n ) should by relatively a 
big since the system is invariant under translations (mod n) and under rotations. So, one may 
expect P(l,k n — 1) ~ 3n 2 . There exists several programs which calculate this polynomial from 
the coefficients of the system of inequalities which define the polytope. One of these programs is 
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called LattE and it written by Matthias Koppe and Jesus A. De Loera. This method remains to 
be implemented and investigated in a different project. Also, it seems like the number of variables 
that can be used in this program is also by about 100. However, we think that the method fits 
very well with the toroidal formulation. If the Ehrhart polynomial can be computed in dimensions 
k and 2k, it may give more information of how to find the best density. We will be using other 
methods which are discussed in the next sections. 

3. Semi-regular tessellations 

Working with the graphs generated by semi-regular tessellations is more challenging when it 
comes to implementing the problem into LPSolve. There are also advantages here since the systems 
have less symmetry and somehow this is a plus for the optimization programs to arrive faster at a 
maximum. 

3.1. The case T = (3 3 ,4 2 ). In Figure [6] (a) we see an arrangement of a half-domination set with 
a density of X, which is only t4 off of the upper bound in the next theorem. 




(a) Density (b) Density 

Figure 6. Some arrangements for T= (3 3 ,4 2 ) 



Theorem 3.1. The half- domination density for the tessellation T = (3 3 ,4 2 ) satisfies 

13 

Proof. Let us consider a toroidal graph induced by G-j^ m n with m and n big and a half domi- 
nation set corresponding to it. We observe that there are = mn squares and |(3mn) = 2mn 
equilateral triangles in Gf^ m ^ n . As usual let us consider the variables x v for each vertex v in this 
graph defined to be 1 or as being in the domination set or not. Also, we denote by x* v the sum 
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of the variables x w corresponding to the adjacent vertices w of v. For a vertex v corresponding to 
a square we have (see Figure [6]) 



(4) 



2x v + x* v <4. 



We denote by T the sum of all x v over all vertices corresponding to triangles and by 5 the sum of 
all x v over all vertices corresponding to squares. If we sum up all equalities @ over all the squares 
we get: 



25 + J2x* v < Amn =► 25 + (25 + T) < 4mn or 45 + T < 4mn. 

v for a squares 

For a vertex v corresponding to an arbitrary triangle we have 



(5) 2x v + x* v <3, 

which gives, as before, if summed up over all of the triangles: 



2T + x* v = 6mn 

v for a trianlge 



2T + (25 + 2T) < 6mn or 25 + 4T < 6mn. 



If we let x = jj^- and y = jp^j, we need to maximize x + y and, as we have seen above, x and y 
are subjected to the restrictions 



4x + y < | 
2x + 4y < 2. 

The usual maximization technique (see Figure [7|) gives that for real numbers x and y, the maximum 
of x + y is attained for x = and V = m with x + y = This proves that 



P(3 3 ,4 2 ) 



lim sup 

m.n— >oo 



S(m, n) + T(m, n) 
3mn 



< 



13 
21' 



3.2. Toroidal examples for T = (3 3 ,4 2 ). The number of variables grows rapidly with m = n. 
For n = 7 we have 3n 2 = 147 variables and this seems to be a good bound for what one can 
obtain with LPSolve. We use variables x^j^ (having values or 1) with i = 0, 1,2, ...,n — 1, 
j = 0, 1, 2, m — 1 and k = 1, 2, 3 with Xj j i for a square, £ij,2 and £ij,3 for the triangles in the 
cluster of minimal tiles as in the Figure [71 The system of inequalities can be written as follows 
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0.2 0.4 0.6 




x(0,0,1) 



x(1,0,1) 




Figure 7. The two lines and the variables convention 



2 x i,j,i + x i,j,2 + x u,j,2 + x v,j,i + x i,t,3 < 4, where u = i — 1 (mod n), v = i + 1 (mod n), 
< w = j — 1 (mod m), £ = j + 1 (mod m), u,v £ {0, 1, 2, n — 1}, w,t £ {0, 1, 2, m — 1} 
k 2 ^i,i,2 + x i,j,l + x u,j,3 + x v,j,3 < 3, and 2xjj i3 + x ijj2 + x U j )2 + < 3. 

The best densities we have gotten are listed next: 
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7 

12 


3 

5 


11 

18 


a* 



Of these arrangements, the case m = n = 6 (Figured] (b)) gives the highest density which is within 
Y^g to the upper bound shown above. 

The same method used before for finding upper bounds does not give a better bound as the one we 
have proved in Theorem 13. 21 For m = n = 7 we let 2^0,2 = ^j,o,3 = and xoj,i = x o,j,2 = x o,j,3 = 
with i € {0, 1, 2, ...,n — 1}, j G {0, 1, 2, m — 1} and LPSolve gives a maximum of 72. So, the 
upper bound is 147 _ (7( j+ 7(3) _ 2) = § > i • 

In [7] we have introduced the concept of deficiency function, Sij k, and global deficiency A of an 
arrangement. Let us see how this works in this situation. We define 
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if ■Ei,j,k 



2 if k = 1 

1 if k = 1 or 2, 



and 



if x 



v -<j,fc — 3 if A; = 1 or 2, 
A = — — ^2 ^j,j,fcj where V is the set of vertices. 

We observe that the arrangement in Figure[S](a) has A = = — -I and the arrangement in FigureE] 
(b) has A = — yg. We point out that the closer the global deficiency, A, is to zero, the bigger the 
density of an arrangement is. It seems like an arrangement in which A = is not possible. 

3.3. The case T = (3, 6, 3, 6). As before we will use the same technique to find the "trivial" upper 
bound which is now usual for these half domination problems. This upper bound is also sharp in 
this case. 





(a) Density 



(b) The infinity column 



Figure 8. "Trivial" arrangement 



(c) Domain 



For a vertex v corresponding to a hexagon we get 3x v + x* < 6, and for v corresponding to a 
triangle we have 2x v +x* < 3. We denote by H the sum of all x v over all vertices corresponding to 
hexagons and by T the sum of all x v over all vertices corresponding to triangles. These inequalities 
imply 



3H + (2H + 2T) < 6mn and 2T + {AH + T) < 6mn. 
The system in x = and y = becomes 

' 5x + 2y < 2 

< 

Ax + 3y<2, 
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and the maximum of x + y is equal to | and attained for x = and y = | (see Figure M (c)). 
Nevertheless, there seems to be a surprising resemblance between this optimal solution and the well 
know sculpture of C. Brincu§i. (see Figure [3] (b)). 



3.4. The case T = (3, 4, 6, 4). The system in this case becomes a little more complicated since we 
have three different type of tilings, but the gives an upper bound which is strictly less than |. 

Theorem 3.2. The half- domination density for the tessellation T = (3, 4, 6, 4) satisfies 

19 

P(3,4,6,4) < 

Proof. The inequalities defining the problem are given by 



2x v + x* < 3 for v corresponding to a triangle 



2x v + x* < 4 for v corresponding to a square 



3x v + x* < 6 for v corresponding to a hexagon. 



As before we introduce T = Y x v , S = Y x v, and H = Y x v- The 

v for triangle v for square v for hexagon 

inequalities above give 



2T + 25 < 3(2mn) 



2S + (3T + 6H) < 4(3mn) or < 



3H + 2S < 6mn, 



2x + 2y < 1 



3x + 2y + 6z < 2 



2y + 3z < 1, 



where x = y = and z = g^j- The usual optimization methods give the maximum for 

x + y + z, under the above constrains and x,y,z € [0, 1), to be attained for x = g, y = ^ and 



z = and a value of 



The variables we are going to use for writing the system for LPSolve are indexed as before JCij fc, 
withi G {0, l,2,...,n-l}, j G {0, 1, 2, m - 1} and k G {1,2,3,4,5,6} (see Figure Ej) . Given 
as above, we let u,v G {0, 1, 2, n — 1}, and w,t E {0, 1, 2, m — 1} such that u = i — 1 (mod n), 
u = i + 1 (mod n), w; = j — 1 (mod m), and t = j + 1 (mod m). Then the system can be written 
in the following way (see Figured]): 
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Figure 9. Variables convention 



(6) 

^•Cjj. I ~\~ %i,j,2 %i,j,4 %i,j,6 %u,jfi ~i~ %i,w,4 ~i~ %v,w,2 ^ 6 
2^Cj,j,3 Xi,j,2 ~l" x u,j,4 ~\~ %u,t,6 — 3, 2X{jfi + 2^i,j',4 %i,j,6 -^v,j,2 — 3, 

2xij j2 + Xij,i + Xi,j,3 + %u,j,5 + Xu,t,i < 4, 2xij i4 + Xij t i + + Xijf + < 4, and 

Contrary to what is expected LPSolve takes more time to solve these systems even under one 
hundred variables. The best densities we have gotten are listed next: 



(m,n) 


(1,1) 


(2,2) 


(3,3) 


(4,4) 


/ 5 (3,4,6,4),(m,n) 
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2 


7 

12 


31 

54 


7 

12 




( a ) P(6,4,3,4) 



( b ) P(6, 3,3,3,3)— 9 

Figure 10. Some distributions 
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3.5. The rest of the semi-regular tessellations. We are going to include for the rest of the 
semi-regular tessellations only the most significant facts found but without proves. One case use the 
same methods to check them. For the semi-regular tessellations (8 2 ,4) and (12,6,4) the following 
arrangements (Figure [TlT a) and (b)) gives the best densities. 




(a) P( 8 2 i4) = | (b) P(12,6,4)=§ 

Figure 11. Best arrangements 



We notice that the deficiency for each of the arrangements in Figure [TlT a) and (b) is equal to 
zero. 




0) P(l22,3) = f (P) p(4, 3,3,4,3)=| 

Figure 12. Best arrangements 



For the semi-regular tessellation (12 2 ,3), we have gotten the theoretical upper bound of | and 
one of the arrangements as in Figure fT2T a) which we think it is actually the sharp. In Figure fT2Tb). 
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we have a sharp arrangement for tessellation (4, 3, 3, 4, 3). Finally, for T = (6, 3, 3, 3, 3) we see see 
an arrangement of density | in Figure [10] (b), but one can show that the best density is actually | 
given by the distribution shown in Figure [5J 

4. Conjectures and other comments 

From what we have seen so far, there are some patterns that emerge. Given a vertex transitive 
infinite graph (for every two vertices, there exists a graph isomorphism mapping one vertex into 
the other) have half-domination arrangements which have rational best densities. It is not clear 
if such arrangements are unique (up to the isomorphisms of the graph) or there exist essentially 
different variations. In any case, we see that if the deficiency is zero, then the solution seems to 
be unique. If the deficiency is positive, one may expect to have more solutions and we have such 
an example in the case of the King's Graph (see [6]). The bigger the deficiency the higher the 
number of combinatorial possibilities that can result in maximum arrangements but, we conjecture 
that there are only finitely many of them. Results that show the exact number of such maximum 
arrangements are, nevertheless, at our interest in further investigations. However, we believe that 
the right methods to approach these questions successfully, even with the assistance of powerful 
computers, are yet to be discovered. Another path of investigations is to look into finding similar 
answers to ^-dependence problems in all of the graphs studied here. 
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